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Abstract 

We present a complete analysis of the electric dipole moment (EDM) of the 



> 

m 

P^ . p-meson induced by CP violating operators of dimension 4 and 5 within the 

^^ . QCD sum rules approach. The set of CP-odd operators includes the theta 

Qv \ term and the electric and chromoelectric dipole moments of the three light 

On \ quarks. We find that the /?- meson EDM induced purely by the EDMs of 

quarks is smaller, but still in reasonable agreement, with the predictions of 

a naive constituent quark model. However, the chromoelectric dipole mo- 

Mh. ments, including that of the strange quark, give comparable and sometimes 

r~| ! larger contributions. We also consider the effect on the hadronic EDM of the 

existence of Peccei-Quinn symmetry. When this symmetry is active, chro- 

L^ ■ moelectric dipole moments induce a linear term in the axion potential which 

j_j ■ leads to a numerically important vacuum contribution to the hadronic EDM. 
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I. INTRODUCTION 

Tests of time-reversal symmetry at low energies are an important source of information 
about the CP properties of the physics at and above the electroweak scale, complementary to 
that coming from K and B meson physics. Impressive experimental progress achieved during 
the last decade has brought the limits on the electric dipole moments (EDMs) of neutrons, 
heavy atoms, and molecules |l|] down to a remarkably low level. The Kobayashi-Maskawa 
model, so successful in explaining the observed CP violation in K mesons, predicts EDMs 
to be several orders of magnitude smaller than the current experimental sensitivity. This 
presents a unique opportunity for limiting extra sources of CP-violation, and the constraints 
resulting from EDM data are generally very strong [|2|. 

Generically, EDMs can be used to probe the physics at a high energy scale by limiting 
the coefficients of operators Oi with dimension A; > 4 in the effective low energy Lagrangian. 
The effective Lagrangian for these operators has the form, 

C.ff -Y.<^M'-'0f\ (1) 

i 

where M is the mass scale at which these effective operators are induced and Cj their co- 
efficients which, in general, have logarithmic scale dependence. These operators are odd 
under CP transformations and their coefficients q are proportional to the fundamental CP- 
violating phases of the underlying theory. Consequently, the calculation of a hadronic EDM 
can naturally be separated into two main parts. Firstly, there is the calculation of the 
coefficients Cj for a specific model of CP-violation which involves integrating out distances 
shorter than M~^. The second part, which is by far the more complicated, is the problem 
of switching from the perturbative quark-gluon description to the level of hadrons which 
requires nonperturbative input. 

In this letter, we present a systematic analysis of the EDM of the simplest hadron - the 
p-meson (more specifically p"*") - induced by operators with dimension 4 and 5 within QCD 
sum rules which up to now remains the most reliable analytical method for calculating the 
properties of hadrons [^. We choose the following parametrization of dimension four and 
five CP violating sources: 

SC = e-^GG - - J2 dq-qFa-i^q - - H dqqgsGa-f^q, (2) 

q=u,d,s q=u,a,s 

where F^^ and G^^ are the electromagnetic and gluonic field strength tensors. The first term 
here is the effective theta term, which primarily is due to the fundamental QCD vacuum 
angle Oqcd- This represents a challenge for particle physics model building as naturally one 
would expect ^ ~ 0{1) which is ruled out by all available data on EDMs. As a consequence, 
one is usually led to introduce the Peccei-Quinn (PQ) mechanism via which this primary 
source of 6 is removed, Oqcd = 0. Nonetheless, even in the presence of PQ symmetry, as 
we shall discuss, the effective 6'-term is non-vanishing. The two sums in (^ represent the 
EDMs and chromoelectric dipole moments (CEDMs) of the three light quark flavors. Other 
quarks are considered to be heavy and can thus be integrated out producing operators of 
dimension 6 and higher. 



To motivate this calculation, we note that in principle experiments [0] impose strong 
constraints on some combination of the coefficients 6, di, and di and thus on the funda- 
mental CP-phases of the theory. In order to extract these constraints, however, we have 
to embark on a non-perturbative calculation of the electric dipole moment of the neutron 
in particular (or nucleon-nucleon interaction, as in the case of the EDM of ^^^Hg) induced 
by the effective Lagrangian (^. The connection between different EDM observables and 
the coefficients 6, di, di is especially important in the framework of supersymmetric theo- 
ries where eq. (^ and an additional three-gluon CP-odd operator [Q represent a complete 
set of the relevant operators with the coefficients explicitly calculable as functions of the 
soft-breaking parameters. (For a recent discussion in the context of the MSSM, see e.g. 0). 

Generically, the EDM of a hadron can be written as a linear combination of the coeffi- 
cients in eq. (^ 

4 = dh{9) + rf^°^(rf., dd, 4) + ci^''°''(rf~., d,, ds). (3) 

In the case of the neutron EDM, dn{0) has been calculated using various different tech- 
niques: (1) making use of the dominance of a pion loop-induced logarithm in the chiral limit 
10]; (2) in the Skyrme model [^; and finally using QCD sum rules PJT(][]; and all have pro- 
duced similar results. For the calculation of d^^^{du, d^, dg), the quark tensor charges over 
the nucleon are required, and various techniques lll|iIH have produced results consistent 
with the predictions of a naive SU(6) quark model. We note that the contribution of the 
strange quark EDM is found to be consistent with zero in all methods. 

Unfortunately, the quantitative evaluation of the CEDM contribution to the neutron 
EDM is considerably more complicated and although a number of serious attempts have 
been made to estimate d'^^^^{du,dd), these results often differ by more than one order of 
magnitude. In particular, direct QCD sum rules calculations |1T3|JT^ give dn ~ 20 times 



smaller than the estimates based on the chiral loop approach |]T3|,|T^. Thus, we believe 
that an independent calculation of d^^^^{du,dci,ds) within the QCD sum rule approach is 
absolutely necessary in order to clarify the magnitude and sign of the contribution of quark 
CEDMs to the EDM of the neutron. 

With this motivation in mind, the discussion of the p-meson presented hereQ will serve 
a dual purpose. Firstly, the detailed analysis of the contribution of dimension four and five 
CP-odd operators to a particular hadronic EDM allows us insight into the relative sizes 
of the contributions, and in particular the dominant mechanism inducing the EDM, which 
may well prove quite universal. Secondly, the p-meson is the simplest light quark system 
where QCD sum rules are known to work well. Thus it is a convenient arena in which to 
develop techniques that should later prove useful in an analogous, but much more involved, 
study of nucleon sum rules and the neutron EDM. In this regard, we recently calculated 
dp{6) by generalizing the operator product expansion (OPE) to the case of an external 6 
background |T^, which was subsequently applied to the case of dn{0) 0. In this letter this 
method is developed further to include all CP-odd dimension 5 operators. As a result we 
obtain dp as an explicit function of 6, di, and rfj. We observe the numerical importance of 



^Note that the p-EDM was also studied within the context of a quark model 



all CEDM contributions and a strong dependence of the result on the existence/absence of 
PQ symmetry. With the use of the techniques developed here and in our earlier work, the 
question of d^^^^{du, d^, dg) will be readdressed in the near future |T9 |. 



II. OPE ANALYSIS 

Within the sum rule approach to the p"*" EDM, we need to consider the correlator of 
currents with p"*" quantum numbers, in a background with a CP violating source and an 
electromagnetic field F^,^, 

n,.(V) =^ J d'xe'^-{0\T{j;{x)j;m\0)qP,F. (4) 

In the presence of these sources, and since we will work outside the chiral limit with un- 
equal quark masses, it is necessary to take into account mixing between the vector current 
associated with p"*" and the axial- vector current [|TB|. Thus we parametrize the full current 
in form 

j; = V^ + tcA^, (5) 

where c is a (real) mixing parameter to be determined, and 

VfM = u'j^d A^ = M7^75(i. (6) 

The parameter c is linear in the CP violating source, and thus to first order we have, 

n,. = {V,Vj)cp + ^c{{V,Al) - {A,Vj)) + ■■■ (7) 

In evaluating these correlators, we require the propagator in the presence of CP violating 
sources, 5*^^ = {qa{x)qi,{0))cp,F, which has the form 

a xe b^^,{x) - 1 — — — — — . [fi) 

Using eq. (P), we can express the relevant current correlators as 

' d\e'^-%Q\V^{x)V}{<dm = -^ f d'xe'P-^TTi^,S''{x)^,S^i-x)), (9) 



with similar expressions for (VA^) and {AV''). The appropriate tensor structure to isolate 
is F^u{= ^fj.uafsF"''^ /2), and we evaluate the correlators as usual by projecting to the vacuum 



in order to extract this structure (see e.g |]T8[ for more details in the case of p"*"). In doing 



so we will assume a constant background electromagnetic field, and use a fixed point gauge 
||20| for the gluon field. Condensates involving quark and gluon fields are then parametrized 
in terms of certain "condensate susceptibilities" following loffe and Smilga |^ , 

gs{0\qGcrq\0) = -ml{0\qq\0) 
{0\qa^.q\0)F = egxF^,{0\qq\0) 
gsiOlqiG^nqlO)^ = e,KF^,{0\qq\0) 
2^.(0|g75(Gp'^)g|0)F = te,^F^,{0\qq\0). (10) 
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FIG. 1. Contributions to the correlator at leading order in Ff^y. 

In general, the OPE will contain perturbative pieces in addition to terms explicitly propor- 
tional to the quark-gluon condensates. However, it turns out that the simple perturbative 
contributions are heavily suppressed. This is due to the chirality-flipping structure of the 
EDM and CEDM operators. The result then vanishes unless there is an additional chirality 
flip in the loop which can only come from a quark mass in the propagator. Indeed, for the 
6'-induced contribution, the perturbative piece is identically zero 



Thus, we are left with only the non-perturbative contributions, for which the leading 
order diagrams are given in Fig. 1. Naively, one could drop all the terms in the quark 
propagator which are proportional to a small mass, m„ or md- This is not valid, however, as it 
turns out that at different stages of the OPE calculation, quark masses from the propagator, 
m„(^), enter in front of vacuum correlators of the form C ~ / d'^x{Q\T{u'~^^u, u{Ga)'~^^u{x)\Q) . 
These correlators can be saturated by pion exchange, producing a non-vanishing contribution 
in the chiral limit, niqC ~ 0{mq/m\) — > const. Equivalently, one can perform a chiral 
rotation, requiring that pions cannot be produced from the vacuum, (0|£e//|7ro) = 0. This 
leads to the appearance of additional 75-mass terms in the Lagrangian which will certainly 
contribute to dp. We prefer, however, not to perform the chiral rotation and instead to 
account for the vacuum contributions explicitly. 

Using this approach, we find for {VV"^), 

(^m(p)K^(O)) = — F y'^ [mdeuiu-f^u) c/P - m^ed{d-/r^d) c/p) 

+ {dd - du) (qq) - ( /« " 9 ) *^^'^^" ^ ^"^'^'* ^^^'^ ) ' ^^^^ 

where the first term arises from Fig. 1(a), the second from Fig. 1(b), and the third from 
Fig. 1(c). A useful simplification in these calculations follows by noting that, although 
naively contributing to the structure F^j^^, via the equations of motion, contributions arising 
from a Taylor expansion of the quark wave function, actually vanish! Thus we find the 
rather compact expression exhibited in Eq. (|ll]). 

The necessity for additional contributions, due to mixing with the axial vector current, 
now becomes clear. Indeed, if eq. (0) were the entire answer, the first line calculated in the 
external 6 background would produce an incorrect quark mass dependence, differing from 



mumd/ {rriu + rud) ■ Thus its clear that we need to address mixing with the axial current |]T8 
The mixing parameter c may be obtained to leading order by considering the correlators 
(VV'^ — AA^), and {VA^ — AV"^) with the external field turned off. At this order we may 
diagonalize on the tensor structure gp.yp^, and we obtain 



-{dd — dy) ; 

6 niu + md 



— i- 






(12) 



This value of c ensures the absence of mixing to leading order between the "eigen-currents" 
Vfj_ + icA^ and iA^ + cV^. The calculation of the mixed correlators is straightforward, and 
the result is given by 



lyA^) = -{AV^) 



^u 



pz 



X(m„ed + mrfe„). 



(13) 



Combining the relevant pieces according to eq. (0) we arrive at the following result: 



n 



tF, 



/lU 



flU 



pZ 



{qq) 



dd- du + xi^d - eu] 



m^md {d'y^d + wy^u) 



CJP 



rriu + md 



Xm^idd - du 



7 N eurrid + edTUu 



rriu + md 



+ K 



{ddCu - duCd) 



(14) 



The next problem to address is the calculation of the vacuum matrix elements in eq. 
(|l^). These terms, arising from Fig. 1(a), require the evaluation of correlators of the form 
J d'^y{q'-y^q{x),i6C{y)), where 6C given in Eq. (|) involves in particular the ^-term and the 
colour EDM sources which may be extracted from the vacuum at leading order in the 
background electromagnetic field. The case of 9 was discussed at length in our previous 
paper ||18| . Here we shall concentrate on the CEDMs. We will evaluate these correlators by 
inserting a complete set of intermediate states {ej}. 



{q-f5q{x),iSC{y)) = {q-i^q{x) \Y.\^i) 2^^^! ^^^iv))- 



—■rrii 



(15) 



For the two-flavour case, to a good approximation we need only consider ttq, while for three 
flavours we also include t]. The neglected states then have masses > IGeV and will remain 
massive in the chiral limit. They may contribute corrections of no more than 25%, which 
we shall factor into our precision estimate. The remaining condensates may be reduced to 
commutators in a manner analogous to the soft pion theorem in chiral perturbation theory. 
For a generic operator O we have 



(0|0|vro) = -V(0|[O,MV-ciSrf]|0) 



2U 



{0\O\v) 



2V3U 
Concentrating on the colour EDM sources in 



{0\[O,u^-f5U + d^j5d-2s^-f5s]\0). 



we then obtain 



(16) 

(17) 



(m75m)^o 

(M75M)r, 



■{d'J^d) 



TTO 



-{du - dd){qq) ^ 

2 mu + nid 



(rf75rf)r, = --:ds{qq)^- + 
4 m» 



(18) 
(19) 



where in the last expression we have neglected terms of 0(c?(„^rf)/ms), assuming an approxi- 
mate proportionahty of the CEDMs to the quark masses, i.e. dd/dg ~ md/ms <^ 1. 

Putting the pieces together, and decomposing into singlet and triplet combinations, we 
finally obtain the following result for the correlator in (0), 



n 



tF, 



^lU /_ 



/iy 



{m) 



'—e_m^,^ + d_ 
2 






e+(i_) 



(20) 



where we have introduced the notation: d± = du ^ dd, d± = du i: dd, e± = e^ ± e^, 
m^ = mumd/{mu + m^), and m = (m„ - md)/{mu + nid). 

III. THE EDM AND PECCEI-QUINN SYMMETRY 

The phenomenological side of the sum rule for the correlator II^i, may be obtained from 
the form-factor Lagrangian, C = —2idpmpp't F^i,p^ + ■ ■ ■. We find. 



n 



Phen 



-2id„m„F, 



A2 



p-" fj-l' 



(j)2 _ rj^'iy 



+ 



(21) 



where A is the coupling to p^, and we have kept only the double pole term. We now 
follow standard practice, equating this result with pO| ) and performing a Borel transform to 
suppress the contribution of continuum states. Rather than presenting a detailed analysis 
of the sum rule, including a parameterization of single pole and continuum states, we shall 
consider only the double pole term here, as previous work p!8|,|9| indicates that the corrections 
which arise from a more careful analysis are overwhelmed by errors from other sources. At 
this order, the coupling may be obtained in terms of the Borel mass M as A^ = m^M^/47r^ 



from the CP even sum rule (see e.g. [^). It turns out that within these approximations dp 
does not depend on M^. We then obtain the following result for the EDM induced by 6, 
{du, dd) and the colour sources {du, dd, dg), 



^EDM+CEDM+e 



27r2 



ni'i 



■{m)) 



+ a_ H — 7::^a_(e+ — e_mj 



X"^o T m, 1 / 
2 m^ 2 \ 



12 



(e^djL. — e^d- 



(22) 



This result parameterizes the effect of all the dimension four and five sources, including the 
^-term which was considered previously in []T8|. 

The numerical values for the condensates were obtained in [^], |Q and |p!^ . 



m. 



0.8 GeV^ [23] 

-5.7 ±0.6 GeV"2 |2| 



-0.34 ±0.1 14 



i = -0.74 ±0.2 [III 



(23) 
(24) 
(25) 
(26) 



Note that with these values the combination (k — ^/2) actually vanishes within the specified 
precision. Substituting the numerical values for the condensates, masses, and charges, we 
obtain a final result for the EDM of the p-meson: 

4 4 . in-3 _ _ _ '? 'iMeV ~ 
dp = '^^^^ e + 0.51(rf, - dd) - 0.13e(d„ - da) - —^ e4 (27) 

Here we have used m„ = 4.5 MeV, and md = 9.5 MeV. We will comment on error estimates 
in Section IV. 

It is interesting to observe that the central point for the QCD sum rule prediction for 
dp induced purely by the EDMs of the quarks is essentially half the prediction of the naive 
non-relativistic quark model, 

d^^QM = d^-dd (28) 

Nonetheless, given the relatively large error bounds on (|27|) , and noting that smaller values 
of mu and irtd will bring the two results closer together, we feel that the results are in 
reasonable agreement. 

The contributions of up- and down-quark CEDMs are not, in fact, any smaller than those 
of the quark EDMs. This may be seen in a simple example of supersymmetric CP-violation 
due to a squark-gluino loop. The induced EDM operators will then be proportional to 
the charges of the quarks, thus bringing in an additional factor of e/3 to the contribution 
induced by dd, leading to a result of the same order as the CEDM contributions. Our result 
for (ip(CEDMs) is larger than one would expect from power-counting rules such as "naive 
dimensional analysis" |2^ which are often applied to estimate (i„ (CEDMs). The importance 
of CEDM contributions in the QCD sum rule approach follows directly from the values of the 
X and especially m^ condensates. It seems clear that the estimate within naive dimensional 
analysis corresponds to a group of perturbative diagrams in which the gluon from dq is 
attached to a quark line. This would correspond to an 0(as/(47r)) suppression relative to 
the leading order OPE terms considered here. 

Perhaps the most interesting point to make here is the significance of the contribution 
arising from the strange quark CEDM operator. The l/irts suppression is in fact fictitious 
as we expect dq ~ m^. When this ansatz is assumed, the contribution of the s-quark CEDM 
becomes even larger than that of the down quark! 

However, the presence of the theta term is expected to be numerically dominant, as it 
corresponds to a dimension four operator and does not experience any suppression by a 
heavy mass scale. This poses a serious problem for any model of CP-violation other than 
Kobayashi-Maskawa, as corrections to the theta term are normally large and need to be 
cancelled by extreme fine-tuning of the "initial condition" for theta. Currently, the accepted 
recipe to avoid a "^-dominated" EDM is to assume Peccei-Quinn symmetry. This mechanism 
is apparently a necessity for any SUSY model operating with CP-SUSY phases around the 
electroweak scale. 

The existence of this mechanism brings an additional contribution to the EDM, not 
contained in eq. (^2]). PQ symmetry, although allowing the axion to set 6'qcd = 0, still 
leads to CP violating terms due to linear contributions to the axion potential pGf. In 
particular, the axion potential has the form V ~ —9^K — 26K', where 



K = i(^GG,^GG) (29) 

is the topological susceptibility, and 

K' = I (^GG, -lrf^gGa75g) (30) 

are correlators arising from the CEDM sources. This linear shift in the axion potential then 
leads to an "induced" ^-terni with coefficient 

,„, . _^ . !| j: Jl. (31) 

q=u,d,s 1 



This result follows by evaluating the correlators in the manner described in [B^] (see also 



28|). The value of ^jnd is, of course, independent of the axion coupling constant, and of 



the particular manner in which PQ symmetry is implemented. Consequently, we find an 
additional vacuum contribution to the EDM in ([lT|) and ([12D of the form, 

{uj5u)e + {di5d)e = iOmdiqq) = -ml i J2 -^] (^l)- (^2) 

\q=u,d,s 1 J 

Including this expression]^ in ( p2D we observe the complete cancelation of the term propor- 
tional to the strange quark CEDM, while this new vacuum source of du and dd is combi- 
natorially less suppressed than the direct contributions appearing in (p2D. The final result 
then has the form 



^EDM+CEDM+PQ ^ ^(_(^g)) 



m? 



a„ H a_ C-L — 6-171) 

12 ^ + ^ 



H -{d^ — d-m)e_ H — Ik (e_(i+ — e+(i_ 



(33) 



Inserting numerical values for the condensates, masses, and charges as before, we find: 

dp = 0.51{d^ - dd) - 0.34e(J„ - dd) - 0.58e((i„ + dd). (34) 

We see that in the presence of PQ symmetry, the contribution from the CEDM sources is 
even more pronounced. 

IV. DISCUSSION 

We have studied the EDM of the p-meson induced by CP-odd operators of dimension 
four and five. We find that, at leading order, QCD sum rules predict that dp induced by the 



^Note that this contribution only affects terms proportional to x^^o ^^ ^^^ contributions propor- 
tional to ^ and K are subleading in the EDM induced by ||l^ 

9 



quark EDMs is within a factor of two of the prediction of the nonrelativistic quark modeL 
Moreover, a conservative estimate of the error, accounting for higher dimension operators, 
neglected intermediate states, and single pole contributions, would suggest that this factor 
of two could easily be accomodated within the precision of our estimates. This precision 
could of course be improved by including next-to-leading order corrections in the OPE, 
and also accounting for the single pole contributions on the phenomenological side. An 
important advantage of the QCD sum rules approach to the calculation of hadronic EDMs 
over other methods is the proportionality of the result to dq{qq) which is normalization-point 
independent as in most of the models we have dq ~ rriq. The uncertainty related with the 
poor knowledge of nig for the light quarks is thus significantly reduced. 

We find that dp induced by the colour EDM operators is actually comparable in magni- 
tude with dp{EDMs). There is no specific mechanism of suppression which can be attributed 
to (ip(CEDMs). This conclusion supports an estimate of dn made in Ref. ^^. We also wish 
to stress the numerical importance of the color EDM operator of the strange quark, whose 
contribution is of the same order as that arising from the up and down CEDM operators. 

Peccei-Quinn symmetry, which removes 6* ~ 1, is an apparent necessity for any model 
with large CP-violating phases other than Kobayashi-Maskawa. When the PQ mechanism 
is active, the axion vacuum experiences a linear shift induced by CEDMs. We find that the 
contribution of the up and down quark CEDMs in this case becomes even larger, whereas 
the contribution of the strange quark CEDM is completely canceled at this order. This 
is an important result which may well prove generic and apply also to the experimentally 
relevant case of the neutron EDM. This important problem will be addressed in a subsequent 
publication |T9[ . 
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